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We develop the single-channel local potential for the KN system, which is applicable to quanti¬ 
tative studies of K bound states in nuclei. Because the high precision measurement of the kaonic 
hydrogen by SIDDHARTA reduces the uncertainty of the KN amplitude below the KN threshold, 
the local potential should be calibrated in a wide energy region. We establish a new method to con¬ 
struct the local potential focusing on the behavior of the scattering amplitude in the complex energy 
plane. Applying this method, we construct the KN potential based on the chiral coupled-channel 
approach with the SIDDHARTA constraint. The wave function from the new potential indicates 
the KN molecular structure of A(1405). 

PACS numbers: 13.75.Jz,14.20.-c,11.30.Rd 


I. INTRODUCTION 

The multi-nucleon systems with an antikaon draw sig¬ 
nificant attention in hadron and nuclear physics. It is 
considered that the strong attraction in the KN chan¬ 
nel leads to various interesting phenomena. The sim¬ 
plest example is the A(1405) resonance as a quasi¬ 
bound state mm- The difficulty of the description of 
A(1405) by the three-quark picture in the constituent 
quark model [3] is compatible with the interpretation as 
the quasi-bound state of the KN system slightly below 
the threshold. In the early days, A(1405) in the KN 
scattering amplitude was analyzed with i^-matrix ap¬ 
proaches and dispersion relations iiHini- Later, coupled- 
channel approaches with chiral SU(3) symmetry [TTHT5] 
have been applied, and confirmed the KN quasi-bound 
picture of A(1405). In addition, a recent lattice QCD 
analysis supports this picture based on the vanishing of 
the strange magnetic form factor m- The KN molecule 
picture indicates that the spacial structure of A(1405) is 
relatively larger than the usual hadronic scale, which is 
shown by several approaches [nHH]. Experimentally, 
the A(1405) signal in the ttE spectrum has been studied 
with various reactions [22]-[26]. 

Another interesting example is the antikaon bound 
states in nuclei, the ^-nuclei [27H29] . Because of the 
strong KN attraction, the ^-nuclei may exhibit the 
qualitatively different structure from the normal nuclei. 
Experimentally, there have been some claims for the 
evidence of the KNN state [3QH33) . Eor instance, J- 
PARC E27 experiment has reported a broad enhance¬ 
ment in the proton coincidence missing mass spectra in 
the (i(7r+,i^+) reaction at 1.69 GeV/c [33]. However, 
we have to note that quantitative results of these ex- 
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periments are not consistent with each other. Eurther- 
more there are experiments which have found no such 
quasi-bound structure [34f[36| . To draw a definite con¬ 
clusion, further studies are needed. Theoretically, the rig¬ 
orous three-body calculations of the KNN system have 
recently been performed [la nasi. All calculations 
have obtained the qualitatively consistent result that the 
KNN system is bound between the KNN and ttYN 
thresholds. However, the quantitative predictions of the 
mass and the width are substantially different from each 
other and are not consistent with the experimentally re¬ 
ported values. In this way, the quantitatively conclusive 
result of the KNN system has not been achieved. 

It is the KN interaction below the threshold that is 
essential for the calculations of ^-nuclear systems. How¬ 
ever, the subthreshold region cannot be directly accessi¬ 
ble by experiments, so we have to extrapolate the scat¬ 
tering amplitude constrained by the experimental data 
above the KN threshold. Previous studies of the K- 
nuclei have suffered from the large uncertainty, mainly 
because the experimental data has not been sufficient to 
constrain the subthreshold amplitude. Recently, the SID¬ 
DHARTA collaboration has measured the precise energy- 
level shift of the kaonic hydrogen [45l [46] . This data is 
related to the K~p scattering length 07], which quanti¬ 
tatively constrains the scattering amplitude at the KN 
threshold. This reduces the uncertainty of the ampli¬ 
tude below the KN threshold significantly [48l|49]. Eor 
a reliable prediction of A(1405) and the ^-nuclei, the 
constraint from the SIDDHARTA data should be taken 
into account. 

The base for the few-body calculations is the two-body 
hadron interaction. Historically, the hadron interaction 
has been constructed phenomenologically. In the case 
of the nuclear force, phenomenological interactions are 
quite successful in reproducing the experimental data 
with the precision of y^/d.o.f 1 (50] E]. Though 
the phenomenological interactions have been successfully 
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applied to various few-body systems, the direct connec¬ 
tion to QCD is not obvious. The first principle calcu¬ 
lation of QCD is the lattice simulation which provides 
the promising approach to the hadron potentials [5^ . 
However, the nuclear force in the realistic set-up is yet 
to be constructed. Another approach is based on chi¬ 
ral perturbation theory which is the effective field theory 
of QCD with chiral symmetry being the guiding prin¬ 
ciple IS3 EH- In this approach, the potential can be 
systematically improved with the higher order contribu¬ 
tions. In the state-of-the-art calculations, it is possible to 
construct the nuclear force as precise as the phenomeno¬ 
logical ones. 

In this work, we construct the KN potential using chi¬ 
ral unitary approach which is based on chiral perturba¬ 
tion theory and unitarity of the scattering amplitude. 
Thanks to the systematic improvement, the low energy 
K~p total cross sections, threshold branching ratios, and 
the SIDDHARTA data are well reproduced with an ac¬ 
curacy of x^/d.o.f ^ 1 [48l [49]. The potential is con¬ 
structed in the local form in the coordinate space for 
the convenience of the applications to few-body calcula¬ 
tions. In contrast to the nuclear force, the KN potential 
cannot be directly obtained in chiral perturbation theory 
which does not contain the long range meson exchange 
processes. We therefore construct the potential so as 
to reproduce the scattering amplitude from chiral uni¬ 
tary approach on the real energy axis following Ref. [55| . 
Given that the uncertainty of the subthreshold amplitude 
is reduced by the SIDDHARTA constraint, we have to 
establish the construction procedure with the high pre¬ 
cision in the wide energy region. Moreover, to analyze 
the structure of A(1405), the precision in the complex 
energy plane is necessary. In this way, we construct the 
reliable KN potential applicable for the quantitative cal¬ 
culations. 

In Sec.jlT] we briefly introduce chiral unitary approach 
for the KN scattering, and the framework to construct 
the hadron local potential from this approach. In Sec. Ing 
we examine the construction procedure to reproduce the 
original amplitude even in the complex energy plane with 
a simple model as an example. The new construction 
procedure of the hadron potential is applied to the KN 
amplitude with the SIDDHARTA constraint in Sec. |IV[ 
leading to the reliable KN potential. Using this new KN 
potential, we investigate the spatial structure of A(1405). 
The last section is devoted to the summary of this work. 


US] which is based on the resummation of the interaction 
terms derived from chiral perturbation theory. The 5- 
wave meson-baryon scattering amplitude Tij{y/s) at the 
total center of mass energy y/s is 
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( 1 ) 


where Uj and Gi represent the meson-baryon interac¬ 
tion kernel derived from chiral perturbation theory and 
the loop function, respectively with the meson-baryon 
channel indices being denoted by i,j. There are four 
meson-baryon channels with isospin / = 0, KN, ttS, 77 A 
and KE corresponding to i = 1,2,3 and 4, respectively. 
The interaction kernel Vij is systematically obtained in 
chiral perturbation theory, where the leading contribu¬ 
tion is given by the Weinberg-Tomozawa term. System¬ 
atic improvement with higher order correction has been 
discussed in Ref. [in m ESHSS]. Recently, the refined 
calculations for the S = —1 sector including the next- 
to-leading order terms [48l |49l l6Qti6^ are being available 
with the constraint from the SIDDHARTA data. The di¬ 
mensional regularization is applied to the loop function 
Gi with the finite part being specified by the subtraction 
constant. Adjusting the subtraction constant adequately, 
the experimental data such as scattering cross sections, 
threshold branching ratios and the scattering length can 
be reproduced well. Although there are other regulariza¬ 
tion schemes constrained by the crossing symmetry m 
and the SU(3) symmetry [63], the present phenomeno¬ 
logical regularization scheme is sufficient to consider the 
KN scattering near the threshold. 

The KN forward scattering amplitude is related 
to the amplitude Tij as 


^kn{V^) — ~ 


4:7r^/s 




( 2 ) 


where represents the nucleon mass. In the isospin 
/ = 0 channel, there are two resonance poles in the 
A(1405) energy region, induced by the attractive inter¬ 
actions of the KN channel and the ttS channel [551164]. 
In this paper, we refer to the higher (lower) energy pole 
near the KN (ttS) threshold as KN pole (ttS pole). 


B. Equivalent single-channel potential 


II. FORMULATION 

A. Chiral SU(3) dynamics for KN scattering 

To describe the KN scattering, it is mandatory to con¬ 
sider the channel coupling with the lower energy ttE state 
and the existence of the A(1405) resonance below the 
threshold. Here we utilize the nonperturbative coupled- 
channel framework called chiral unitary approach [TI1 - 


Our aim is to construct the KN single-channel interac¬ 
tion for the application to few-body calculations as well 
as the A(1405) analysis. In this work, we construct a 
single-channel local potential which produces the ampli¬ 
tude equivalent to the chiral coupled-channel approach. 
The coordinate space wave function calculated by the po¬ 
tential is useful to study the spatial structure of A(1405). 
In addition, the local potential is easily implemented in 
the variational calculations of the few-body systems [65] . 
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To this end, we first extract the single-channel KN 
interaction from the coupled-channel scattering equa¬ 
tion 0 . We define the effective interaction as 


N 


N 


^11 = E ^^ruGmVml + E , 


m=2 




md=2 


(3) 


N 


'' ml 


kl 


k=2 


_ q(n-i) 


1 -1 


, m, / = 2, 3, N. 


The quantities with the superscript {N — 1) are the 
{N — 1) X (N — 1) matrices. Using this single-channel 

scattering equation Tn = — Gi] , the origi¬ 

nal amplitude is exactly reproduced. Because of the elim¬ 
ination of the lower energy ttH channel, the effective KN 
interaction V^i has an imaginary part. 

Next, we define the energy dependent local potential 


U{r, E) = g{r)N{E)V^f{E + Mn + tuk), 


N(E) 


Mn (^k + E]sf 

2{E -h Mn + tuk) ujkEn ’ 


(4) 

(5) 


where E, En and ujk are the nonrelativistic energy, the 
energy of the nucleon and the energy of the anti-kaon. 


E = y/s — Mn — ttik^ 
s - 


En = 


luk = 


2v/^ 

s - M% + m\ 


the parameter b by the matching of the full amplitude at 
the KN threshold. This prescription is along the same 
line with the determination of the flux factor N{E). 

The potential well reproduces the original ampli¬ 
tude near the KN threshold, while the deviation in¬ 
creases in the energy region far below the threshold. To 
enlarge the applicability of the potential, we add the cor¬ 
rection AU {E) to the strength of the potential, 

U{r, E) = g{r)N{E) [V,f{E + Mn + m^) + AV{E)] . 

(7) 


For the analytic continuation of the amplitude in the 
complex energy plane, it is useful to parameterize the 
strength of the potential by a polynomial in the energy. 


U{r = 0,E)=g{r = 0)N{E) 



E 


100 MeV 


(8) 


We refer to the energy range where the potential is pa¬ 
rameterized as parameterized range, which will be speci¬ 
fied for each potential. We comment on the analytic be¬ 
havior of the amplitude calculated from the potential ^ . 
Because the potential is constructed to reproduce the 
original amplitude, the correct analytic behavior is guar¬ 
anteed within the parametrized range on the real axis. 
On the other hand, the extrapolation of this potential 
to other energy regions should be carefully performed, 
since some unphysical singularities can in general be de¬ 
veloped. This will be discussed in detail in the next sec¬ 
tion. 


III. POTENTIAL CONSTRUCTION 


with the mass of the antikaon m^. The spatial distribu¬ 
tion of the potential is governed by g{r) which is normal¬ 
ized as J dr g{r) = 1. The flux factor N{E) is determined 
by the matching with the original amplitude at the KN 
threshold in the Born approximation m- In this work, 
we choose a Gaussian for g{r) 


9{r) = 


7r3/2^3 




where the parameter b determines the range of the po¬ 
tential. Using the local potential, we can calculate the 
wave function from the Schrodinger equation. 


In this section, we study how the original amplitude is 
reproduced by the KN local potential. Examining the 
previous method in Ref. m in detail, we improve the 
construction procedure to reproduce the original ampli¬ 
tude even in the complex energy plane. Here, we mainly 
employ the amplitude of the HNJH model Eaini] for the 
comparison with Ref. [55] . Inclusion of the SIDDHARTA 
constraint will be discussed in the next section to con¬ 
struct a realistic KN potential. 


A. Precision of potential in the complex plane 


1 d‘^ ij ( 

= Eu{r), (6) 

where g = Mnttik/{Mn + tuk) is the reduced mass 
and u{r) is the 5-wave part of the two-body radial wave 
function. From the behavior of the wave function at 
r ^ oc, the scattering amplitude E^n can be obtained. 
In Ref. [55], the parameter b was determined to match 
the amplitude Frn with the original amplitude in the 
A(1405) resonance region. In this work, we determine 


A resonance state is represented by a pole of the scat¬ 
tering amplitude in the complex energy plane. The pole 
structure of the KN amplitude is therefore important for 
the study of the spacial structure of A(1405). It is con¬ 
sidered that the pole structure of the KN system may 
affect the result of the KNN system m- We thus focus 
on the scattering amplitude from the previous potential 
in the cormlex plane. 

In Fig. [h we compare the KN (/ = 0) scattering am¬ 
plitude from the local potential E^n R^cf. [55] with 
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FIG. 1. (Color online) Scattering amplitudes from the lo¬ 
cal potentials Frn (thick lines) and the amplitudes directly 
from chiral unitary approach (thin lines) with models 

ORB [68], HNJH BNW [S6l[57] and BMN [58]. The 

real (imaginary) parts are shown by the solid (dotted) lines. 


Model 

pole position [MeV] 

^ KN 

Frn 

ORB [SH] 

1427 - 17z, 1389 - 64z 

1419 - 42z 

HNJH EiEZl 

1428 - 17z, 1400 - 76z 

1421 - 35z 

BNW [571 [59] 

1434 - 18i, 1388 - 49i 

1404 - 46z 

BMN [SD 

1421 - 20i, 1440 - 76z 

1416 - 27z 


TABLE I. Pole positions of the original scattering amplitudes 
from chiral unitary approach and the amplitudes from 

the local potentials Frn- All poles are found in the ttE 
unphysical and KN physical Riemann sheet. The pole at 
1440 — 76i in the BMN model is above the KN threshold, 
and hence is not in the most adjacent sheet to the real axis. 


the corresponding original amplitude in the chiral unitary 
approach F^^ for the models ORB [68], HNJH [661(67] . 
BNW [56] [57] and BMN [58] on the real axis. The KN 
amplitudes on the real axis are reasonably well repro¬ 
duced by the potentials in Ref. [55]. On the other hand, 
we find a large deviation of the amplitude in the complex 
energy plane. In Table [l] we list the pole positions of the 
scattering amplitudes. While chiral unitary approaches 
generate two poles in the A(1405) energy region, the local 
potentials give only one pole. In addition, the position of 
the pole does not agree with neither of the original poles. 
Hence, the potential construction procedure should be 
improved by paying attention to the amplitude in the 
complex energy plane. 

To improve the construction procedure, we introduce 
several quantities to assess the deviation of the am¬ 
plitudes in the complex plane. For the discussion of 
A(1405), we consider that the following energy region 


is relevant!^ 


1332 MeV < Re[z] < 1450 MeV 
— 100 MeV < Im[ 2 ;] < 50 MeV, (9) 

where z represents the complex energy of the two-body 
system. First, we define the average deviation AF^eai 
between the amplitude from the local potential Frn and 
the amplitude from chiral unitary approach F^^ on the 
real energy axis as 

f dy/s |F^^(v^) - FRR{y/s)\ 

AFreal = J -• (10) 

J dyfs \F^\{R8)\ 

When AFreai is small, the amplitude on the real axis 
is well reproduced by the potential. When AFreai ^ 1, 
it means the average deviation on the real energy axis 
reaches the same amount as the average magnitude of 
\^kn\‘ With the HNJH model, we obtain AF^eai = 0.14. 

Second, we define the deviation of the amplitude at 
complex energy 2 ;, 


AF{z) 



( 11 ) 


In this paper, we regard that the amplitude is well repro¬ 
duced when the deviation is smaller than 20 %: 


AF{z) < 0.2. 


( 12 ) 


We call the energy region satisfying this condition the 
“precise region”. We also define the percentage of this 
precise region in the relevant energy region ©by 


comp 


Jj d{Rez)d{lmz) 0(0.2 - AF{z)) 

JJ (i(Re.2)d(Im2:) 


X 100. 


(13) 


If the local potential well reproduces the original ampli¬ 
tude well in the relevant region of the complex energy 
plane, then we have Pcomp ^ 100. The HNJH model 
gives Pcomp = 50, which quantifies the insufficiency of 
the precision in the complex energy plane. 


B. Region near the real axis 

We explain how to reproduce the amplitude in the com¬ 
plex energy plane. We first focus on the region near the 
real energy axis including the KN pole. Here we use 


^ The lower boundary of Re[ 2 ;] (1332 MeV) is set at the ttE thresh¬ 
old. 
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1340 1360 1380 1400 1420 1440 


Re z |MeVl 

FIG. 2. (Color online) The contour plot of AF of the HNJH 
potential in Ref. [55] . The unfilled region corresponds to large 
deviation, AF > 2. Precise region is defined as AF < 0.2. 
The crosses represent the original pole positions of A(1405). 



1340 1360 1380 1400 1420 1440 


Rez IMeVl 

FIG. 3. (Color online) The contour plot of AF of Potential I. 
The unfilled region corresponds to large deviation, AF > 2. 
Precise region is defined as AF <0.2. The crosses represent 
the original pole positions of A(1405). 


the HNJH model as an exampl^ Let us show the con¬ 
tour plot of AF in the complex energy plane with the 
potential in Ref. [55| in Fig. Here we choose the most 
adjacent Riemann sheet to the real energy axis. It is 
seen from Fig. j^that the deviation in the region around 
Re[z] ^ 1400 MeV is larger than the other region. The 
deviation of the amplitude should influence the pole po¬ 
sitions of A(1405). The reason for the deviation is that 
the correction to the potential AV has been applied only 
in the region below 1400 MeV in the previous work. Fur¬ 
thermore, the AV has been chosen to be real, based on 
the dominance of the real part in Vii . 

In this work, we add AV in the relevant energy re¬ 
gion for the KN pole, 1332-1450 MeV. Hereafter we call 
the region where AV is applied the correction range. To 
reproduce the original amplitude near the A(1405) res¬ 
onance region,we introduce the complex correction AV. 
As a consequence, AF^eai is significantly reduced. We 
call the new potential with the complex AV “Potential 
I” and summarize its properties in Table [HI together with 
the property of the corresponding potential in Ref. [55] , 

With Potential I, the deviation on the real energy axis 
AFreai IS rcduccd by two orders of magnitude. Thanks to 
the reduction of AFreai, the KN pole position is also sig¬ 
nificantly improved. We show the contour plot of AF in 
Fig|^ Comparing Fig. Hwith Fig.|3l we find that the pre¬ 
cise region {AF < 0.2) of Potential I satisfying Eq. (12) 
is extended over the KN pole. The improvement of the 
pole position can be understood by this enlargement of 


the precise region. Quantitatively, Pcomp in Eq. (13) in¬ 
creases from 50 to 68. In this way, the precision near the 
real axis can be improved by introducing the complex 
correction AV in the relevant correction range. 


C. Region far from the real axis 

While Potential I reproduces the original amplitude 
near the real energy axis, the deviation of the amplitude 
increases in the region far from the real axis (see Fig. 
and the ttE pole does not appear. Here we further im¬ 
prove the potential paying attention to the region far 
from the real axis. 

In principle, if the original amplitude is completely re¬ 
produced in the whole range on the real energy axis, 
the analytic continuation in the complex energy plane 
is unique. This suggests that the increase of the param¬ 
eterized range will inmrove the precision of the potential 
far from the real axisjj On the other hand, there is a lim¬ 
itation of extension of the parameterized range because 
of the threshold effect. In the present framework of the 
effective single-channel potential with polynomial param¬ 
eterization, it is difficult to incorporate the non-analytic 
threshold effect of the other channels. The parameterized 
range can only be extended to the nearest thresholds. 
In this case, the parameterization of the KN potential 
strength should be performed between the ttE threshold 
(1331 MeV) and the r]A threshold (1664 MeV). In order 


^ In this work, the range parameter b of the potential is deter¬ 
mined as 0.46 fm by the new prescription explained in Sec. ^ ^ In this subsection, the correction range is chosen to be the same 

in contrast to b =0.47 fm of the potential in Ref. m with the parameterized range. 
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Ref. [55] 

Potential I 

Potential II 

b [fm] 

0.47 

0.46 

0.46 

AV 

real 

complex 

complex 

polynomial type 

third order 

third order 

tenth order 

correction range 

1300-1400 

1332-1450 

1332-1521 

parameterized range 

1300-1450 

1332-1450 

1332-1521 

AFreal 

1.4 X 10“^ 

4.8 X 10“^ 

4.0 X 10“^ 

Pcomp 

50 

68 

85 

pole positions [MeV] 

1421 - 35i 

1427 - 17i 

1428 - 17i, 1400 - 77i 


TABLE IL Properties of the HNJH potential in Ref. [55] and Potential I and Potential II in this work. Shown are the potential 
range parameters 5, the corrections to the strength of the potentials AE, the polynomial types of the potential strength in 
energy, the correction ranges where AV is applied, the parameterized ranges by the polynomials, the average deviations AEreai 
from the amplitudes of chiral unitary approach on the real energy axis, the percentages of the precise region in the complex 
energy plane, and the pole positions of the amplitudes from the potentials Frn- The pole positions of the original amplitude 
are 1428 - 17i MeV and 1400 - 76i MeV. 


to keep the precision on the real axis for the larger param¬ 
eterized range, we increase the degree of the polynomial 
from the third order to the tenth order. 

To examine the above strategy, we construct the po¬ 
tentials varying the parameterized range by 1 MeV. The 
typical results of AFreai, ^comp, and the pole positions 
of these potentials are shown in Table |ml In all cases, 
AFreai IS rcduccd by an order of magnitude from that of 
Potential I. This is because we change the parameteriza¬ 
tion from the third order to the tenth order polynomial. 
Though the wider fitting range leads to the slightly larger 
AFreai, the Order of magnitude remains same. In general, 
when a high-degree polynomial is used for the parame¬ 
terization, artificial poles appear between the KN and 
ttE thresholds. In the present case, this occurs when the 
fitting range is smaller than ^ 1500 MeV. However, as 
the fitting range increases, these unphysical poles move 
away from the relevant energy region and only two phys¬ 
ical poles remain. The KN pole appears at the original 
pole position, 1428 — 17i MeV and is stable against the 
variation of the parameterized range. On the other hand, 
the position of the ttE pole depends on the parametrized 
range. The optimized value of the upper boundary of the 
parameterized range is 1521 MeV to reproduce the orig¬ 
inal pole position, 1400 — 76i MeV. At the same time, 
the maximum value of Pcomp is achieved. We call the 
potential with the best parametrized range Potential II. 
We show the contour plot of AF with Potential II in 
FigJ^ As shown in Figj^ we succeed in extending the 
precise region to Im^; ^ —80 MeV, near the ttF pole. As 
a consequence, we obtain two poles both at the correct 
positions. 

It turns out that the largest parameterized range does 
not always lead to the best potential. In the present case, 
this is because the ttF pole position moves along with the 
change of the parameterized range. The best potential is 
achieved when the ttF pole comes closest to the original 
position. 



1340 1360 1380 1400 1420 1440 


Re 2 IMeVl 

FIG. 4. (Color online) The contour plot of AF of Potential II. 
The unfilled region corresponds to large deviation, AF > 2. 
Precise region is defined as AF <0.2. The crosses represent 
the original pole positions of A(1405). 


IV. APPLICATION 


In the previous section, we have established the con¬ 
struction procedure to reproduce the original amplitude 
in the complex energy plane, considering the high preci¬ 
sion on the real energy axis and the wider parameterized 
range. In this section, we apply this procedure to chiral 
unitary approach with SIDDHARTA constraint [481149] 
and construct the realistic KN local potential. This new 
potential is then used to estimate the mean distance be¬ 
tween K and nucleon, that is, the spatial structure of 
A(1405). 
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upper boundary [MeV] 

APreal 

P 

-L comp 

pole positions [MeV] 

1450 


59 

1428 - 

17i, 1388 — 60z, unphysical poles 

1500 

2.6 X 10“"^ 

71 

1428 - 

17i, 1404 — 70z, unphysical poles 

1521 

4.0 X 10“"^ 

85 

1428 - 

17z, 1400- 77i 

1550 

5.9 X 10“^ 

79 

1428 - 

17z, 1392 - 82i 

1600 

6.8 X 10“^ 

77 

1428 - 

17i, 1389 - 83z 

1650 

8.8 X 10“"^ 

77 

1428 - 

17i, 1389 - 87i 


TABLE III. Results of the precision of the potentials against the variation of the parameterized range for the HNJH model. 
Shown are the average deviations AEreai from the amplitudes of chiral unitary approach on the real energy axis, the 

percentages of the precise region in the complex energy plane, and the pole positions of the amplitudes from the potentials 
Frn- The “unphysical poles” stand for the artificial poles generated between the KN and ttE thresholds as explained in the 
text. The pole positions of the original amplitude F^^ are 1428 — 17i MeV and 1400 — 76z MeV. 


A. Realistic KN potential 


As we explained in Sec. [I| the constraint from the pre¬ 
cise SIDDHARTA data is crucial for the quantitative cal¬ 
culation of the K and nucleons systems. In this section, 
we construct the KN local potential based on the ampli¬ 
tude of Refs. [IHl |49] with the SIDDHARTA constraint. 
To apply to the few-body A-nuclei, we construct the po¬ 
tential of the 1 = 1 amplitude in addition to the / = 0 
channel. 

The amplitude of Refs. [48l |49] is given in the parti¬ 
cle basis with the isospin breaking effect in the hadron 
masses. On the other hand, the potential in the isospin 
basis with isospin symmetry is useful for various applica¬ 
tions. Moreover, in the practical potential construction 
procedure, the existence of multiple thresholds in the par¬ 
ticle basis prevents us from enlarging the parameterized 
range. We thus construct the isospin symmetric KN am¬ 
plitude by replacing the physical hadron masses by the 
isospin averaged ones keeping the low energy constants 
and subtraction constants the same as Refs. [48l|49]. The 
result of the isospin symmetric KN amplitude (/ = 0) 
is shown in Fig. [^together with the combination of the 
original amplitudes 

of Refs. [48l [49]. The difference stems from the isospin 
breaking effect. From this figure, we find that these am¬ 
plitudes well agree with each other except for the tiny 
region near the KN threshold. Since the difference in 
the most important region for A(1405) and the KNN 
systems is negligible, we adopt this isospin symmetric 
amplitude to construct the KN potential. 

Following the construction procedure in Sec. |n| here 
we determine the Gaussian parameter b = 0.38 fm. We 
show the properties of the potentials with various param¬ 
eterized ranges in Table |TV| The optimal upper boundary 
of the parameterized range to reproduce the pole posi¬ 
tions is found to be 1657 MeV. We call the best potential 
SIDDHARTA potential (7 = 0). The properties of the 
SIDDHARTA potential (/ = 0) are summarized in Ta¬ 
ble |V| and the contour plot of AF is shown in Fig. 
We find that SIDDHARTA potential (/ = 0) well re¬ 
produces the original amplitude in the complex energy 



FIG. 5. (Color online) The result of the isospin symmetric 
KN amplitude in the 7 = 0 channel (thick lines) and the com¬ 
bination of the original amplitudes (TV-p^-p + 27V-pxo^ + 
Fk^uK^ti) with isospin breaking [48l[l9] (thin lines). The 
real parts are shown by the solid lines, and the imaginary 
parts are shown by the dotted lines. 


plane (AF^eai = 5.4 x 10“^, Pcomp = 96), and the poles 
of A(1405) appear at the same position in the accuracy of 
1 MeVThe strength of the potential is shown in Fig. 
as a function of the energy. The energy dependence of the 
potential strength is not strong, but is important to pre¬ 
cisely reproduce the original amplitude. The coefficients 
of the strength Ki in Eq. are shown in Table 

In the same way, we construct the KN local poten¬ 
tial for the 1 = 1 channel from the combination of 

The range param- 

eter of the potential is determined to be 5 = 0.37 fm. In 
this channel, however, the natural analytic continuation 
of the amplitude is not possible, because of the prescrip- 


VI 


^ We note that the maximum of Pcomp is achieved when the upper 
boundary is set to be 1658 MeV. Because the value of Pcomp 
depend on the definition of the relevant region (|^, we determine 
the best potential by the accuracy of the pole positions. 












upper boundary [MeV] 

AEreal 

P 

-L comp 

pole positions [MeV] 

1450 

0.91x10"^ 

50 

1424 — 28i, 1381 — 49z, unphysical poles 

1500 

1.7x10"® 

62 

1424 — 26z, 1395 — 62i, unphysical poles 

1550 

2.4x10"® 

70 

1424 - 26z, 1379 - 68i 

1600 

2.8x10"® 

72 

1424 - 26z, 1381 - 70z 

1650 

4.8x10"® 

86 

1424 - 26z, 1382 - 79z 

1657 

5.4x10"® 

96 

1424 - 26z, 1381-81i 


TABLE IV. Results of the precision of the potentials against the variation of the parameterized range for the amplitude with 
SIDDHARTA constraint. Shown are the average deviations AEreai from the amplitudes of chiral unitary approach on the 
real energy axis, the percentages of the precise region in the complex energy plane, and the pole positions of the amplitudes 
from the potentials Frn- The “unphysical poles” stand for the artificial poles generated between the KN and ttE thresholds 
as explained in the text. The pole positions of the original amplitude F^^ are 1424 — 26z MeV and 1381 — 81i MeV. 


SIDDHARTA potential (/ = 0) 

b [fm] 

0.38 

AV 

complex 

polynomial type 

tenth order 

parameterized range [MeV] 

1332-1657 

AEreal 

5.4 X 10“^ 

Ecomp 

96 

pole positions [MeV] 

1424 - 26z, 1381 - Sli 


TABLE V. Properties of SIDDHARTA potential (/ = 0). 
Shown are the potential range parameters b, the corrections 
to the strength of the potential AV, the polynomial type of 
the potential strength in energy, the parameterized range by 
the polynomials, the average deviation AEreai from the am¬ 
plitude of chiral unitary approach F^^ on the real energy 
axis, the percentage of the precise region in the complex en¬ 
ergy plane, and the pole positions of the amplitude from the 
potential Frn- The pole positions of the original amplitude 
are 1424 - 26i MeV and 1381 - 81z MeV. 


tion to avoid the unphysical cut of the amplitude [57] . In 
contrast to the / = 0 channel, the best value of the upper 
bound of the parameterized range cannot be determined 
from the information of the complex energy plane. Here 
we use the same parameterized range as that in the / = 0 
channel. This may be sufficient for the present purpose 
because the interaction in this channel is not as strong 
as the / = 0 channel and the contribution to few-body 
systems is considered to be small. In this way, we con¬ 
struct SIDDHARTA potential (/ = 1) whose strength at 
r = 0 and the coefficients Ki are shown in Fig. and 
Table |Vn| respectively. As expected, the strength of the 
real part of the potential is smaller than the / = 0 coun¬ 
terpart. The imaginary part is similar magnitude with 
/ = 0, suggesting the absorption occurs equally in / = 0 
and 1=1. 


B. Spacial structure of A(1405) 

We have succeeded in constructing the new KN local 
potential reliable for the quantitative calculations of the 



1340 1360 1380 1400 1420 1440 


Re 2 IMeVl 

FIG. 6. (Color online) The contour plot of AF of SID¬ 
DHARTA potential (/ = 0). Precise region is defined as 
AF < 0.2. The crosses represent the original pole positions 
of A(1405). 


K nuclei. In this section, as a direct application of this 
new potential, we estimate the KN distance to under¬ 
stand the spacial structure of A(1405). 

Generally, the unstable states are expressed as the 
poles of the scattering amplitude in the complex energy 
plane. As an analogy of a bound state, the spacial struc¬ 
ture of an unstable state is reflected in the wave func¬ 
tion at the pole energy. With the solution of the radial 
Schrodinger equation Uz{r) at the complex energy the 
KN wave function in s wave is written as 


^z{r) = 


1 u^{r) 


(14) 


The wave function of a resonance state diverges at r ^ 
oc. Hence the wave function cannot be normalized by 
the standard normalization, 



( 15 ) 
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Kq [fm] Ki [fm] K 2 [fm] [fm] 

Re Im Re Im Re Im Re Im 

-10.833 -1.8149 -2.7962 -0.64315 -0.47980 0.88991 -0.64480 -0.55225 

Ka [fm] i ^5 [fm] Kq [fm] K 7 [fm] 

Re Im Re Im Re Im Re Im 

0.44645 0.0048399 0.089658 0.47326 -0.23222 -0.38284 0.027650 -0.072843 

Ks [fm] i ^9 [fm] Kiq [fm] 

Re Im Re Im Re Im 

0.059123 0.22152 -0.024071 -0.099375 0.0022208 0.014415 


TABLE VI. Coefficients Ki in Eq. ^ of the strength of SIDDHARTA potential (7 = 0). 


Ko [fm] Ki [fm] K 2 [fm] 

Re Im Re Im Re Im 

-6.2261 -1.7382 -2.1909 -0.63300 -0.37668 -0.082052 

K 4 [fm] K 5 [fm] Kq [fm] 

Re Im Re Im Re Im 

2.9791 0.52183 -0.53283 0.29745 -3.0760 -0.59389 

Kg [fm] Kq [fm] Kio [fm] 

Re Im Re Im Re Im 

0.64668 0.36800 -0.59141 -0.17606 0.10746 0.026128 


Ks [fm] 

Re Im 

-0.14782 -0.18206 

K 7 [fm] 

Re Im 

1.5430 -0.062035 


TABLE VII. Coefficients Ki in Eq. ^ of the strength of SIDDHARTA potential (7 = 1). 



y/s [MeV] 


EIG. 7. (Color online) Strength of SIDDHARTA potential 
(7 = 0) U{r, E) at r = 0. The real part is shown by the solid 
line, and the imaginary part is shown by the dotted line. 



EIG. 8. (Color online) Strength of SIDDHARTA potential 
(7 = 1) U{r, E) at r = 0. The real part is shown by the solid 
line, and the imaginary part is shown by the dotted line. 


Alternatively, the wave function of the non-Hermitian 
problem can be normalized with the Gamow vector la¬ 
beled by the index G [69l |70] , 



= 1 - ( 16 ) 


nary part, the wave function converges at r ^ 00 Hence 
both the prescriptions (15) and (16) are applicable (see 
also Appendix A). 

As explained in Appendix B, for a problem with an 
energy-dependent potential, we should modify the nor¬ 
malization condition to ensure the conservation of the 
norm and the orthogonality relation between two states. 


In the present problem, the poles of A(1405) are in the 

physical Riemann sheet of the KN channel. Because the ^ In the coupled-channel formulation, the wave function of the ttE 
corresponding eigenmomentum has the positive imagi- channel diverges at r ^ 00 . 
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The modified normalization condition is 
dU{r,E) 


/ 


dr [1 — 


dE 


= 1 . 


(17) 


The expectation value of an operator should be modified 
in the similar way. For comparison, we calculate both 
the expectation values with Eq. (16) and Eq. ( pT| ), and 
label the latter one by the index “mod”. 

The result of the wave function normalized with 
Eq. (16) is shown in Eig. Here we use the precise KN 
pole energy, 1423.97 — 26.28i MeV, in order to achieve 
the enough convergence at r ^ 10 fm. We note that the 
wave function has an imaginary part where the phase is 
uniquely determined by the normalization (16). With 


this wave function, we calculate the expectation value of 


as 


(?’^)g = y r'^'tpGirf. (18) 

The result of the mean squared distance (the root mean 
squared distance) of the antikaon and the nucleon is 
{r‘^)G = 0.79 - 1.2H fm, = 1-06 - 0.57i fm). 

Similarly, the distance with the modified normalization 
condition ([TtI can be calculated as 


^^2) mod = J dr 


1 - 


dU{r, E) 

as 




(19) 


The result of the mean squared distance (the root 
mean squared distance) is = 0.71 — 

1.26i fm, (y^= 1.04 —0.61i fm). It turns out that 
the modification of the normalization condition does not 
change the quantitative result very much. 

In Table VIII[ we compare these results with the pre¬ 
vious estimations. In Ref. m, the radius is calcu¬ 
lated by the form factor of A(1405) in the chiral unitary 
model HU IS]. The result in Ref. [21] is obtained by the 
KN wave function in the complex scaling method with 
the coupled-channel potential model. In both cases, the 
leading order Weinberg-Tomozawa interaction is used, 
without the constraint by the SIDDHARTA data. The 
present result from the next-to-leading order chiral inter¬ 
action with the SIDDHARTA constraint quantitatively 
confirms the results of the previous works. 

If the spacial extent of the A(1405) wave function is 
sufficiently larger than the potential range, the radius 
is determined mainly by the tail of the wave function. 
Because the tail is related to the eigenenergy, we can 
estimate the spacial extent from the eigenenergy E (see 
Appendix A) 




^/2k? 2^J—(iE 


= 0.85 - 0.58i fm. (20) 


In this case, the value is same order as the from 

the wave function. This means that the KN distance is 
sufficiently larger than the range of the potential. 



FIG. 9. (Color online) I — ^ KN wave function ipc with 
the normalization ( [T^ at the A(1405) pole energy, 1423.97 — 
26.28z MeV. The real part is shown by the solid line, and the 
imaginary part is shown by the dotted line. 


y/{r^)G [fm] 

SIDDHARTA potential 

1.06 - 0.57Z 

SIDDHARTA potential (modified) 

1.04-0.61i 

Ref. [20] 

1.22 - 0.63i 

Ref. [21] 

1.22 - 0.47i 

estimation from eigenmomentum 

0.85 - 0.58z 


TABLE VHI. Averaged KN distances at the A(1405) energy. 
For comparison, we show the estimations from the form fac¬ 
tor of A(1405) [20] and from the coupled channel potential 
model [22. 


Though (r^)G' or give us some information 

about the spacial structure, it is not straightforward 
to interpret the complex number. As explained in Ap¬ 
pendix A, the dumpling of the wave function outside the 
potential range is related to the standard expectation 
value with normalization (15), 


(r^) = J dr r^|'0(r)p. (21) 


Here we regard this quantity as the measure of the KN 
distance. As explained in Appendix B, the modification 
of the norm due to the energy dependence of the poten¬ 
tial cannot be applied without using the Gamow vector. 
Therefore, we calculate the KN distance with Eq. (21). 


This is convincing, because the values of ^et^Jg and 


VHI 


The 


^ ( 7 ’ 2 ^mod almost Same as shown in Table 
result of the KN distance with Eq. (21) is found to be 


= 1.44 fm. 


( 22 ) 


Considering the charge radii of the proton and K~ are 
about 0.85 fm and 0.55 fm [71], we find that the KN 
distance is relatively large in comparison with the usual 
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FIG. 10. (Color online) KN density distribution p in Eq. ( [23t 
(thick solid line), the real part of SIDDHARTA potential (/ = 
0) (thin solid line) and the imaginary part (dotted line) at the 
A(1405) pole energy, 1423.97 — 26.28z MeV. 


b [fm] 


VN) [fm] 

0.2 

0.96 - 0.58z 

1.35 

0.38 

1.06 - 0.57Z 

1.44 

0.4 

1.07- 0.57i 

1.48 

0.6 

1.18-0.57Z 

1.57 

0.8 

1.29 - 0.57Z 

1.67 


TABLE IX. The results of the a verage KN distance with the 
Gamow vector normalization ^{r‘^)G and with the standard 
normalization method against the change of the po¬ 

tential range b. 


hadron size. Therefore we conclude that A(1405) is the 
molecular state of the antikaon and the nucleon. In order 
to visualize the spacial structure, we define the density 
distribution 


p{r) = r2|t/)(r)| , 


(23) 


which is shown in Fig. The substantial distribution 
exists outside the potential range 6 = 0.38 fm. 

Finally, we investigate the KN distance against the 
change of the potential range b. The strength of the 
potential is adjusted to rep rodu ce the o rigin al amplitude 
for each b. The results of {r ‘^)q and \/(r^ are shown in 
In all cases, the values of \/(r^ remain larger 
typical hadron size. We find the qualitative 


IX 


Table 
than tHe 


picture of the molecular state is valid irrespective of the 
potential range. 


tion procedure, we have paid attention to the two steps: 
the precision in the complex energy plane and the con¬ 
straint from the recent SIDDHARTA data. This new 
potential is useful for the quantitative calculation of the 
interesting systems such as K few-body systems and the 
A(1405). 

We first establish the procedure of potential construc¬ 
tion by improving the previous work [55]. The previ¬ 
ous potential almost reproduced the original amplitude 
on the real energy axis, while we have found that there 
is a substantial deviation in the complex energy plane, 
including the poles of A(1405). Therefore we need to 
improve the potential construction procedure to repro¬ 
duce the original amplitude even in the complex energy 
plane. We find that the reduction of the deviation on the 
real energy axis AFreai in the wide parametrized range 
is important, based on the uniqueness of the analytic 
continuation in the complex energy plane. Thanks to 
these improvements, we have succeeded in reproducing 
the original amplitude in the drastically large region in 
the complex plane including the two poles of A(1405). 

Next, we have applied the new procedure to the ampli¬ 
tude with the SIDDHARTA constraint to construct the 
realistic KN potential. Here we produce the KN am¬ 
plitude with isospin symmetry from the coupled-channel 
chiral model in Refs. [48l|49]. Based on these amplitudes, 
we have constructed the realistic KN local potentials for 
the / = 0 and 1 = 1 channels. The / = 0 local poten¬ 
tial reproduces the original amplitude in the complex en¬ 
ergy plane including the poles of A(1405). At the present 
time, this is the most reliable local KN potential for the 
quantitative calculations. Applying this new potential 
to A(1405), we have estimated the spacial structure of 
A(1405). The mean distance of K and N is found to be 
1.44 fm. This result shows the meson-baryon molecular 
nature of A(1405). 

As a future perspective, the calculation of the KNN 
system with the new reliable potential is of particular im¬ 
portance. We hope that this result will bring new insight 
in the theoretical and experimental studies of the KNN 
system. 
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V. SUMMARY 


APPENDIX A : AVERAGE DISTANCE OF 
QUASI-BOUND STATE 


We have constructed the new KN local potential (SID¬ 
DHARTA potential) which reproduces the scattering am¬ 
plitude from the chiral SU(3) dynamics. In the construc- 


In this Appendix, we consider the mean squared radius 
of a quasi-bound state in comparison with a bound state 
in the zero range limit. When the spatial extent of the 
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wave function is much larger than the potential range, we 
can treat the potential in the zero range approximation. 
In this case, the mean squared radius (r^) is determined 
only by the eigenmomentum k, where k = in {k > 0) 
for the bound state and k = j > 0) for the 

quasi-bound state. 

In the zero range limit, the radial wave function of the 
(quasi-)bound state in s wave is written as 


u{r) ^ 


(24) 


where u{r) is related to the wave function as tp{r) = 
u(r)/(V47rr). The normalization condition determines 
the factor A. We consider the following two normaliza¬ 
tion conditions. 



1 dr |t/)(r)|^ = 1, 

(25) 

{^P\^P)g = 

1 dr 'ip{r)‘^ = 1, 

(26) 


where the former is the standard normalization whereas 
the latter uses the Gamow vector labeled by G. In the 
zero range limit, these conditions are expressed as 


C 

m) ^ / 

Jo 


j -2Im[fe]r _ 1^1 _ 1 

“ 2Im[fc] “ 




These integrals converge for the bound state or the quasi¬ 
bound state because lm[k] > 0 The normalized wave 
functions are written as 


tp{r) 


A 

R 


lm[k] _ -0 
27r r 


Im[/c] 

27r r 




(27) 

(28) 


where 0 is an arbitrary real constant. In the standard 
normalization, physical observables are independent of 
the phase of the wave function, so 0 is an irrelevant phase. 
In the case of the bound state k = iK,^ Eq. (27) and 


Eq. (28) are equivalent. On the other hand, these wave 


functions of the quasi-bound state are in general different, 
V’ 7^ i’G- 

With these wave functions, we can calculate the mean 


The resonance wave function (Im[fc]< 0) can be normalized only 
by Eq. ( |26| > with appropriate prescription m- 


squared radius. 


(r^) = J dr 


—An r dr 
27r io 

1 


2{lm[k]y ’ 

{r‘^)G = y r'^'tpcirf 


(29) 


—ih /•“ 
—An / dr 
27r Jo 


2{—iky 


(30) 


In the case of the bound state k = in {k > 0), both the 
normalizations give the same result, 

(^ 2 ) ^ ^ ^ (bound state). (31) 


Eor the quasi-bound state, Eq. (29) and Eq. (30) are 
different: 


(r^) = — 

' ' 2k2’ 

{r^)G = 


1 


2 k ,‘^ + dmy — 27^ 


(32) 

(quasibound state). (33) 


Eq. (32) gives a real number for the mean squared ra¬ 
dius, while Eq. (33) gives a complex number. It is com¬ 
mon to use the Gamow vector pi[7ni[72] for unstable 
states. Because the radial wave function of the quasi¬ 
bound state asymptotically behaves as Eq. ( [3Q| ) is 
the natural extension of the bound state. We therefore 
use the normalization with the Gamow vector in Eig. 
and Table lyllTl 

On the other hand, it is difficult to extract the spacial 
information from the complex {r^)G in Eq- (33). We 
note that the dumping of the wave function of the quasi¬ 
bound state is expressed by in the asymptotic 

behavior In this sense, we consider that the real (r^) 
with the standard normalization, which is determined by 
Im[/c], can be interpreted as the spacial extent of the 
quasi-bound state. Hence, in this paper, we use (r^) to 
estimate the spacial extent of the KN quasi-bound state, 
A(I405). 


APPENDIX B : ENERGY DEPENDENT 
COMPLEX POTENTIAL 

As explained in Refs. [TSHIS], the careful treatment 
is necessary for the system with the energy-dependent 
potential. Here, we explain the treatment in the cases of 
real potential and complex one. 

Eirst, we summarize the case of the energy-dependent 
real potential, following Ref. We consider the 
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Schrodinger equation with a time-dependent wave func¬ 
tion ^(r,t)|^with /i = 1, 




at 


-1^^ + Vir4) 


(34) 


For an eigenfunction of the Hamiltonian, = 


becomes 


e the time-independent Schrodinger equation 


H^E{r) = 


lv^ + V{r,E) 


ipE{r) = EipE{r). (35) 


With Eq. (34), the continuity equation for energy- 
dependent potential can be calculated as follows, 


—P = 

at 


at 


4'e + 


a^E 

at 


+V{r,E')^^ 

+ ^E' 




-ii^-^V^ + V{r,E)^^E 

= -V ■ i + i'^%, [V{r, E') - V{r, E)] <i>E, (36) 

where 

P = ^*E,{r,t)'^E{r,t), 

3 = [n'{r,t)V^E{r,t) - {V4/^,(r,f)}v|/£:(r,f)] . 

For the energy-independent potential, the second term of 
the last line in Eq. ( [^ disappears and the usual conti¬ 
nuity equation dPjdt = —V • j can be hold. However, 
for the energy-dependent potential, the additional term 
have to be included. Using the Schrodinger equation, 
id'^Eldt = we obtain the relation. 




V{E') - V{E) 




d^*E, 

at 


E' -E 
V{E') - V{E) 






E' -E 
'V{E')-V{E) 


= {-iE'^E'Y 


E' -E 


V{E') - V{E) 


a^E 

at 




E' -E 
V{E') - V{E) 
E' -E 


4' 


E 


{-iE'^e} 


= i<i!*E, [V{E') - V{E)] ^E, 


(37) 


and the continuity equation for the energy-dependent po¬ 
tential can be modified as 


|(P + P,) = -V^i. 


(38) 


^ Here, we assume that the wave function is normalizable. 


where 

Pa = -n'{r,t) 


V{r,E')-V{r,E) 


E' -E 


4'£;(r,f). 


Therefore, taking the limit of E' 
be modified as 

aV{r,E) 


N = J dr ^*E{r,t) 


1 - 


= J dr -ipEir) 


1 - 


aE 

dV{r,E) 


E, the norm N can 

^Eir,t) 


ipE{r). 


(39) 


Furthermore, the orthogonality relation can be modified 
as 

nr.E')-Vir.Ey 


1 - 


E' -E 


'Ipsir) = 0, (40) 

(E' ^ E). 

Actually, the usual orthogonality relation is not satisfied 
because the term with Pa remains nonzero after integrat¬ 


ing Eq. (38) with respect to r. 

The above method cannot be directly applied to the 
case of the complex potential. In this case, following the 
same procedure. Pa is obtained as 

V%r,E') -V{r,Ey 


Pa = -^yyr,t) 


E' -E 




Because V*{E) y V{E) for the complex potential, tak¬ 
ing the limit of E' ^ E, this term diverges and does 
not become the derivative of the potential. When we 
treat the complex potential problem, it is common to 
use the adjoint wave function T’*' [69], which satisfies the 
Schrodinger equation 




(41) 


Under the adequate boundary condition [69], the time- 

4 


independent eigenfunction -0^ have the following proper¬ 


ties, 

= + ^Yr,E)^ 4(0 = ^*4(0, 


(42) 

(43) 


4(0 = V’£(0> 

^pr,f) = e-*^**4(0 = e“'^*V^(0- 

With this adjoint wave function, we consider the con¬ 
tinuity equation again, labeling the quantities with the 
adjoint function by the index “G”. 

ft* 

E 


9 r.r d^E' . 

_pG ^ - 

at at ^ ^ at 


E 


+ V*{E')\'^ 


4-; 


+ 4* 


-ii--y^ + ViE)}'^E 


= -V • E + (ViE') - V{E)) <i>E, (44) 
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where 


pG 


-I h%{r,t)V^E{r,t) - {V^'l5*(r,t)}^'s(r,i) 


Following the same procedure as in Eq. ( [^ , the continu¬ 
ity equation for the complex energy-dependent potential 
can be satisfied. 


Therefore, we should modify the norm and the orthogo¬ 
nality relation as 


N = J dr dV{r,E) 

= y dr it sir) 


1 

1 - 


dE 

dV{r,E) 


^E{r,t) 

isir), (46) 


^(P® + Pf) = -V-jy (45) 


where 


Pa = -^Ur,t) 


V{r,E')-V{r,E) 
E' -E 




/ 


dr ^E'{r) 


V{r,E')-V{r,E) 
E' -E 


isir) = 0, (47) 
{E' ^ E). 
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